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PROLIFERATING PARASITES IN DIVIDING CELLS: KIMMEL'S 
BRANCHING MODEL REVISITED 

By Vincent Bansaye 

Universite Paris 6 

We consider a branching model introduced by Kimmel for cell 
division with parasite infection. Cells contain proliferating parasites 
which are shared randomly between the two daughter cells when 
they divide. We determine the probability that the organism recov- 
ers, meaning that the asymptotic proportion of contaminated cells 
vanishes. We study the tree of contaminated cells, give the asymp- 
totic number of contaminated cells and the asymptotic proportions 
of contaminated cells with a given number of parasites. This depends 
on domains inherited from the behavior of branching processes in 
random environment (BPRE) and given by the bivariate value of the 
means of parasite offsprings. In one of these domains, the conver- 
gence of proportions holds in probability, the limit is deterministic 
and given by the Yaglom quasistationary distribution. Moreover, we 
get an interpretation of the limit of the Q-process as the size-biased 
quasistationary distribution. 

1. Introduction. We consider the following model for cell division with 
parasite infection. Unless otherwise specified, we start with a single cell 
infected with a single parasite. At each generation, each parasite multiplies 
independently, each cell divides into two daughter cells and the offspring 
of each parasite is shared independently into the two daughter cells. It is 
convenient to distinguish a first daughter cell called and a second one called 
1 and to write + Z^ the number of offspring of a parasite, Z^ of which 
go into the first daughter cell and Z^ of which into the second one. The 
symmetric sharing is the case when (Z( ),Z«) = (zWjZf )). Even in that 
case, the sharing of parasites can be unequal [e.g., when F(Z^Z^ = 0) = 1]. 
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Fig. 1. Multiplication of parasites and cell division. 



We denote by T the binary genealogical tree of the cell population, by 
G n , (resp. G* ) the set of cells at generation n (resp. the set of contaminated 
cells at generation n) and by Z\ the number of parasites of cell i € T, that 
is, 

G„:={0,ir, G::={ieG„:Zi>0}, T := (J G n . 

For every cell i € T, conditionally on Z\ = x, the numbers of parasites 
(Zio, Z\i) of its two daughter cells is given by 

£(zf(i),z«(i)), 
fc=i 

where (zj^\i), (i))} g T,jfc>i is an i.i.d. sequence distributed as {Z^\ Z^) 
(see Figure 1). 

This is a discrete version of the model introduced by Kimmel in [15]. In 
particular, it contains the following model with binomial repartition of par- 
asites. Let Z be a random variable in N and p £ [0, 1] . At each generation, 
every parasite multiplies independently with the same reproduction law Z. 
When the cells divides, every parasite chooses independently the first daugh- 
ter cell with probability p (and the second one with probability 1 — p). It 
contains also the case when every parasite gives birth to a random cluster 
of parasites of size Z which goes to the first cell with probability p (and to 
the second one with probability 1 — p). 

We introduce for a£ {0, 1} 

(1) m a := E(ZW) Vs > 0, f a (s) := E(s zia) ). 
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We assume < tuq < oo, < mi < oo and to avoid trivial cases, we require 

(2) P((Z(°),ZW) = (1, 1)) < 1, P((Z<°\ Z«) G {(1,0), (0, 1)}) < 1. 

This model is a Markov chain indexed by a tree. This subject has been 
studied in the literature (see e.g., [6, 8]) in the symmetric independent case. 
In this case, for every (i, k) G T x N, we have 

P((Z i0 , Z n ) = (k , h) \Zi = k) = F(Z i0 = k \Zi = k)F(Z i0 = k 1 \Z i = k) 

which require that Z<°> and are i.i.d. in this model. Guyon [14] studies 
a Markov chain indexed by a binary tree where asymmetry and dependence 
are allowed and limit theorems are proved. But the case where his results 
apply is degenerate (this is the case momi < 1 and the limit of the number 
of parasites in a random cell line is zero). Moreover, adapting his arguments 
for the theorems stated here appears to be cumbersome (see the remark in 
Section 5.2 for details). In the same vein, we refer to [10, 20] (cellular aging). 

The total population of parasites at generation n, which we denote by 
Z n , is a Bienayme Galton- Watson process (BGW) with reproduction law 
+ Z^K We call Ext (resp. Ext c ) the event extinction of the parasites 
(resp. nonextinction of the parasites), 

Z n = J]] Z{, 

i€<G„ 

(3) Ext = {3nGN:Z n = 0}, 

Ext c = {VnGN:Z n >0}. 

Another process that appears naturally is the number of parasites in a ran- 
dom cell line. More precisely, let (cii)ieN De an i-i-d. sequence independent 
of (Zi)i g T such that 

(4) P(ai = 0) =P(ai = 1) = 1/2. 

Then (Z n ) ng pj = {Z^ ai CL2 v .. jan ))neN is a Branching Process in Random Envi- 
ronment (BPRE). 

The first question we answer here arose from observations made by de 
Paepe, Paul and Taddei at TaMaRa's Laboratory (Hopital Necker, Paris). 
They have infected the bacteria E. coli with a parasite (lysogen bacterio- 
phage M13). A fluorescent marker allows them to see the level of contam- 
ination of cells. They observed that a very contaminated cell often gives 
birth to a very contaminated cell which dies fast and to a much less con- 
taminated cell whose descendance may survive. So cells tend to share their 
parasites unequally when they divide so that there are lots of healthy cells. 
This is a little surprising since one could think that cells share equally all 
their biological content (including parasites). In Section 3, we prove that if 
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D b = {{mo, mi) : morai > 1} 



D\ - {[hi,;. '; : m + mi < 1} 



D 2 - {(m ,mi) : mo + mj - 1} 



D4 = {(mo, mi) : raomi < 1, 

molog(mo) 4- milog(mi) > 0} 



= {(mo. mi) :mo + mi > 1, 
molog(mo) + milog(mi) < 0} 




Fig. 2. 



momi < 1, the organism recovers a.s. (meaning that the number of infected 
cells becomes negligible compared to the number of cells when n — > 00). 
Otherwise the organism recovers iff parasites die out (and the probability is 
less than 1). 

In Section 4, we consider the tree of contaminated cells. We denote by dT 
the boundary of the cell tree T and by dT* the infinite lines of contaminated 
cells, that is 



We shall prove that the contaminated cells are not concentrated in a cell 
line. Note that if mo + mi > 1, conditionally on Ext c , dT* ^ since at each 
generation, one can choose a daughter cell whose parasite descendance does 
not become extinct. 

The rest of the work is devoted to the convergence of the number of con- 
taminated cells in generation n and the convergence of proportions of con- 
taminated cells with a given number of parasites (Section 5). These asymp- 
totics depend on (mo, mi) and we distinguish five different cases which come 
from the behavior of the BGW process Z n and the BPRE Z n (Section 2), 
shown in Figure 2. 

If (mo, mi) E D5, the contaminated cells become largely infected (The- 
orem 5.1). The main two results correspond to cases (mo, mi) E D3 and 
(mo, mi) E D\ and are given by the following two theorems. 

Theorem 1.1. If (mo, mi) E D3, conditionally on Ext c , the following 
convergence holds in probability for every k E N, 



ax = {0,1} 



N 



dT* = {i E dT : Vn E N, Z K \ n ± 0}. 



#{iEG;:Zi = fe}/#G; 



* n— kx> 
n * 



P(T = Jfe) 
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where T is the Yaglom quasistationary distribution of the BPRE (Z n ) n& ^ 
(see [5, 12]). Note that the limit is deterministic and depends solely on the 
marginal laws of {Z^\ Z^) (see Proposition 2.2). This gives then a way to 
compute T as a deterministic limit, although it is defined by conditioning on 
a vanishing event. Kimmel [15] considers the symmetric case ((Z^°\ Z^) = 
(Z^\ Z^)) with mo — m i < 1 < mo + mi in a continuous analogue of this 
model (cells divide after an exponential time) . The counterpart of his result 
in the discrete case is easy to prove [see (21)] and makes a first link with T. 

lim E(#{i £ G n : Z K = k}) /E(#G* n ) = P(T = k). 

Moreover, the proportions of contaminated cells on the boundary of the tree 
whose ancestors at generation n have a given number of parasites converge 
to the size-biased distribution of T letting 77 — > 00 (Corollary 5.4). This gives 
a pathwise interpretation that the limit of the Q-process associated to Z n 
(see [1, 5]) is the size-biased quasistationary distribution. 

Theorem 1.2. // (mo, mi) € D\, (#{i € G* :Z- X = k}) keN conditioned 
on Z n > converges in distribution as n — > 00 to a finite random sequence 

We obtain a similar result in the case (mo, mi) £ D2 (Theorem 5.5) and 
we get the following asymptotics (Theorem 3.1 and Corollaries 5.3, 5.6, 5.8). 

If (mo, mi) 6 D3 (resp. -D5), then conditionally on Ext c , #G*/(mo-(-mi) n 
(resp. #G*/2 n ) converges in probability to a finite positive r.v. 

If (mo, mi) £ D\ (resp. D2), then #G* (resp. #G*/n) conditioned by 
#G* > converges in distribution to a finite positive r.v. 

In the case (mo, mi) 6 -D4, we get only some estimates of the asymptotic 
of #G* which are different from those which hold in the other domains. Our 
conjecture is that #G* has also a deterministic asymptotic, which depends 
on three subdomains (the interior of 7D4 and its boundaries). As a perspec- 
tive, we are also interested in determining which types of convergences hold 
in D4 for the proportions of contaminated cells with a given number of 
parasites (see Section 5.5). 

Moreover, we wonder if the convergences stated above hold a.s. and if 
they extend to the continuous case and complement the results of Kimmel. 
Finally, in a work in progress with Beresticky and Lambert, we aim at de- 
termining the localizations of contaminated cells and the presence of cells 
filled-in by parasites on the boundary of the tree (branching measure and 
multifractal analysis) . 
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2. Preliminaries. In this section, we give some useful results about the 
two processes introduced above. First define: 

(5) m:=\(mo + mi). 

We use the classical notation, where for every i = (ai, . . . , a n ) € G n , 
|i|=n, i\k = (ai, . . . , a^) for every k < n, 
j < i if 3k < n : i\k = j. 

2.1. Results on the BGW process (2 n ) n gj*j. The results stated hereafter 
are well known and can be found in [5]. First, the probability of extinction 
of the parasites satisfies 

P(Ext) = inf{s € [0, 1] :E{s z( ° )+zW ) = s}; 

P(Ext) = 1 iff m + mi < 1/2. 

From now, we assume 

m := E((Z(°> + Z^) log + (Z (0) + < oo. 

Then there exists a random variable W such that 

(6) 7 ^— n -^W, P(W = 0)=P(Ext), E(W) = 1. 

(mo +mi) n 

In the case mo + m\ < 1, there exists b > such that P(i? n > 0) 
b(m,Q + mi) n . Then, there exists U > such that 

(7) P(Z„>0)>C/(m + mir. 

Moreover (2 n ) ng N conditioned to be nonzero converges to a variable called 
the Yaglom quasistationary distribution and we set 

(8) B(s):= lim Ws Zn \Z n >0). 

n—*oo 

We consider then B n k(s) := K(s Zn \ Z n +k > 0) which satisfies 

(9) n lun o ^.( S )= r -^ M . 

Moreover S is differentiable at 1 (Lemma 1 on page 44 in [5]) and we get 

(10) lim \mx B n>k {s) = — U. 

fc^oon^oo <D (1) 

This is the probability generating function of the size-biased Yaglom qua- 
sistationary distribution, which is also the stationary distribution of the 
Q-process. 

Finally if m := E((Z<°) + ZW)((Z(°) + ZW) - 1)) < oo and 2m / 1, then 

(2m) n - 1 



(11) E(Z n (Z n -l))=m(2mY 



(2m) 2 — 2m 
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2.2. Properties of the BPRE (Z n ) n <=^. Recall that (Z n ) n£ ^ is the popu- 
lation of parasites in a uniform random cell line. Then (Z n ) n ^ is a BPRE 
with two equiprobable environments. More precisely, for each n G N, condi- 
tionally ona n = a with a G {0, 1} [see (4)], all parasites behave independently 
of one another and each of them gives birth to Z^ a ' children. The size of the 
population at generation is denoted by k and we note Pfc the associated 
probability. Unless otherwise mentioned, the initial state is equal to 1. For 

the general theory, see, for example, [9, 12, 13, 19]. In the case 
(Z n ) nS N is a BGW with reproduction law Z^°\ 
For i = (ai, . . . , a n ) G G n , we define 

n 

fi ■= fax ° • ■ • fan, "l\ II m <*i> 

i=l 

and for all (n, k) G N x N* and i G G n , 

E k (s z -\(a 1 ,...,a n ) = i) = f i (s) k . 
Then for all (n, k) G N x N* and s G [0, 1], 

(12) E k (s z -) = 2- n f^) k - 

ie<G„ 

First, for every n G N, E(Z n+ i | Z n ) = mZ n and K(Z n ) = m n . 

Moreover, as (P(Z n = 0)) n£ ^ is an increasing sequence, it converges to 
the probability of extinction p of the process. Recalling (1), we have the 
following result (see [19] or [3]). 

Proposition 2.1. Ifmomi < 1, then p = l. Otherwise p<l. 

In the subcritical case (momi < 1), the process Z n conditioned to be 
nonzero which is denoted by Z* converges weakly (Theorem 1.1 in [12]). By 
analogy with BGW, we call its limit distribution the Yaglom quasistationary 
distribution and denote it by T. That is, 

VsG [0,1] E(s Zn |Z„>0)™E(s T ) = C7(s). 

In the subcritical case, the asymptotics of (P(Z n > 0)) nG N when n is large 
depends on the sign of molog(mo) +m 1 log(m 1 ) (see [12]). Now, we require 
also that 

(13) m log(m ) +milog(mi) < 0; E(Z a log + (Z a )) < oo. 

Then we say that Z n is strongly subcritical and there exists c > such that 
as n tends to oo (Theorem 1.1 in [12]), 

(14) F(Z n > 0) ~ cm 11 . 
Moreover, in that case, T is characterized by 
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Proposition 2.2. G is the unique probability generating function which 
satisfies 

G(0) = 0, G'(l)<oo, 

(15) G{h{s)) + G{h{s)) r ,,, n , 
= mG(s) + (1 — m). 

To prove the uniqueness, we need and prove below the following result. 
Lemma 2.3. // H : [0, 1] i-> R is continuous, H(l) = and 

2m 

then H = 0. 

Proof of Proposition 2.2. The finiteness of G'{1) = E(T) is the 
second part of Theorem 1.1 in [12]. 

We characterize now the probability generating function G of the limit 
distribution: 



1 - E(s Zn+1 | Z n+1 > 0) 

_1-E(s z ™+ 1 )_ 1 



■J2nZn = i)(l-Ms Zl )) 



P(Z„ +1 >0) P(Z n+1 >0)^" 

nz n >o) i ( fo( s y+h(sy \ 

¥{z n+1 >^)nzn>Q)h \ 2 / 

F(Z n > 0) 



(1 - E(/o( S ) z " | Z n > 0) + 1 - E(A( S )^ | Z n > 0))/2. 

And (14) ensures that F(Z n > 0)/F(Z n+1 > 0) "-^ m" 1 , so that 
i c() _ l-GWs)) + l-G{A{s)) 

Finally we prove the uniqueness of solutions of this equation. Let G and F 
be two probability generating functions which are solutions of (15). Choose 
a > such that G'(l) = aF'(l). Putting H := G — aF, H' satisfies equation 
(16) and H'{\) = 0. Thus Lemma 2.3 gives H' = 0. As H(0) =0, H = 0. 
Moreover, F(l) = G(l) = 1, so a = 1 and F = G. □ 



PROOF of Lemma 2.3. If H ^ then there exists a G [0, 1[ such that 
/3:=sup{|iI(s)|:sG [0,a]}/0. 



PROLIFERATING PARASITES IN DIVIDING CELLS 9 

Let a n G [0, 1[ such that a n — -h? 1 and a < a n < 1. Then, for every n G N, 
there exists f3 n G [0, a n ] such that: 

sup{|F( S )|: S G[0 ) a„]} = | J H"(/3 n )| 

< |g(/o(/? ra ))|/^(/?n) + \H{h(f3 n ))\f[(f5 n ) 
2m 

<sup{|fl-(a)|:a€[0,l]}, 

since sup{|fl-(s)| : s G [0, 1]} / and (2m)" 1 + /{(/?„))< 1. As /n J = 

0, sup / < sup(I U J) =>■ sup / < sup J, we get 

(3 < sup{|#(s)| : s G [0,a n ]} < sup{|#(s)| : s e]a n , 1]}. 
And H(s) — > leads to a contradiction letting n — ► oo. So if = 0. □ 

In the subcritical case (momi < 1), if molog(mo) + milog(mi) > [resp. 
molog(mo) + milog(mi) = 0], we say that Z n is weakly subcritical (resp. 
intermediate subcritical) and we have P(Z n > 0) ~ c'n~ 3 / 2 -f n [resp. ¥(Z n > 
0) ~ c"n~ l l 2 m n ] for some 7 < m,c' > 0, c" > (see [12] for details). 

Finally we have the following expected result in the supercritical case [4]. 

Proposition 2.4. Ifm m 1 >l, P(Z n *i=3? oo | Vn GN:Z n > 0) = 1. 

3. Probability of recovery. We say that the organism recovers if the num- 
ber of contaminated cells becomes negligible compared to the number of cells 
when n — > oo. We determine here the probability of this event. Actually if 
this probability is not equal to 1, then the parasites must die out for the 
organism to recover. 

Theorem 3.1. There exists a random variable L G [0, 1] such that 

#G*j2 n ™ L. 

If m mi < 1 then P(L = 0) = 1 . 

Otherwise P(L = 0) < 1 and {L = 0} = Ext. 

Remark. In the case mo + mi > 1 and m-omi < 1, the population of 
parasites may explode although the organism recovers. 

This theorem states how unequal the sharing of parasites must be for the 
organism to recover. More precisely, let mo = aM, mi = (1 — a)M where 
M > is the parasite growth rate. Then the organism recovers a.s. iff 

M<2 or a^](l- y / l-4/M 2 )/2,(l + y / l-4/M 2 )/2[ (M>2). 
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Note that for all n G N, 

y(^) = E( ^%: 1Zi>0) =nz n >o). 

Recalling that p is the probability of extinction of (^ n )nGN; 

(17) VneN E^^)=P(Z n >0) n -^l-p. 

The last equality gives also the asymptotic of E(#G* ) as n — > oo in the case 
momi < 1 [see Section 2.2 for the asymptotic of P(Z n > 0), which depends on 
the sign of molog(mo) + mi log(mi)] and in the case m^mx = 1 (see [2, 16]). 

Proof of Theorem 3.1. As #G*/2 n decreases as n increases, it con- 
verges as n — > oo. 

Monotone convergence of ^tG*/2 n to L as n — > oo and (17) ensure that 
E(L) = 1 — p. Using Proposition 2.1, we get P(L = 0) = 1 iff mQ-mi < 1. 

Obviously {L = 0} D Ext. Denote by V(n) the set of parasites at gen- 
eration n and for every p € V(n), denote by Nk(p) the number of cells at 
generation n + k which contain at least a parasite whose ancestor is p. Then, 
for every n £ N, 

{i =o}= n **■<>}. 

P eV(n) 1 ' 

As T n := inf {k > : > n} is a stopping time with respect to the natural 
filtration of (Z{)\\\< n , strong Markov property gives 

P(L = 0) < P(T n < oo)P(L = 0)" + F(T n = oo). 

If P(L = 0) < 1, letting n — > oo gives 

P(L = 0) < n lun o P(T n = oo) = P(Z n is bounded) = P(Ext) 

since i? n is a BGW. This completes the proof. One can also use a coupling 
argument: the number of contaminated cells starting with one single cell 
with n parasites is less than the number of contaminated cells starting from 
n cells with one single parasite. □ 

4. Tree of contaminated cells. Here, we prove that contaminated cells 
are not concentrated in a cell line. If niQ + mi < 1, contaminated cells die 
out but conditionally on the survival of parasites at generation n, the number 
of leaves of the tree of contaminated cells tends to oo as n — > oo. The proof of 
this result will also ensure that, if rriQ + mi > 1, the number of contaminated 
cells tends to oo provided that they do not die out. 
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Theorem 4.1. If mo + mi < 1, #{i € T:Zi ^ 0,Z i0 = 0,Z a = 0} condi- 
tioned by > converges in probability as n — > oo to oo. 
If mo + m-i > 1, conditionally on Ext c , #G* ri — ^ oo a.s. 

Remark. In the conditions of the theorem, (resp. the number of 

leaves) grows at least linearly with respect to n (see Section 5 for further 
results). In the case mo + m\ < 1, conditionally on #G* > 0, the tree of 
contaminated cells is a spine with finite subtrees, as for BGW conditioned 
to survive (see [11, 17]). 

We need two lemmas for the proof. First we prove that the ancestor 
of a contaminated cell has given birth to two contaminated cells with a 
probability bounded from below. We have to distinguish the case where 
P(Z(°) Z^> = 0) = 1, since in that case a cell must contain at least two par- 
asites so that it can give birth to two contaminated cells. 

Lemma 4.2. There exists a > such that for all N G N, i G Gn, n< N 
and fe > 2, 

P(Z j0 ^ 0, Zji + | Zj = fe, Zi > 0) > a 

denoting j = i | n. 

//p(^(o)^(i) = 0) / 1, this result also holds fork = l. 

Proof. We consider first the case ¥(Z^Z^ =0)^1 and we choose 
(Jfco, fci) G N* 2 such that F((Z^°\ ZW) = (fc , fci)) > 0. For every k G N*, we 
have 

P(Z j0 ^ O.Zji ^ | Zj = fc, Zi > 0) > P(Z j0 ^ 0, Zj! ^ | Zj = 1, Zi > 0). 

Moreover, as the function 3mh decreases, we have for all y, x > 

and p G [0, 1[, 

(18) T^>^^- 

1 — maxjy, x j 

Let aG {0,1} and k such that i=jak. Then for all (k ,k[) G N 2 — (0,0), 

P(Z j o = fc ,Z j i = fci |Zj = l,Zi>0) 
P(Zj = ^,Zj 1 = A:i |Zj = l,Z;>0) 

= P(Z(°)=A; ,Z( 1 )=A; 1 |Z ak >0) 

p(z(Q) = ^,za) = ^ |Z ak >0) 

_ P(Z ak > | ZW = kg, = fci) P(ZW = fc , Z^ = fci) 
~ P(Z ak > I Z(°) = fc&, ZC 1 ) = fci) P(Z(°) = fe(„ W) = k[) 
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1 - P(Z k = 0) k " P((Z(°), Z^>) = (k , h)) 
~ 1 - P(Z k = 0) k 'a P((Z(°), ZW) = (A£, jfej)) 

min{fco,fci} P((Z(°),Z( 1 )) = (A; ,fc 1 )) . 
-jfe + fe 1 + ^ + A/ 1 P((Z(o) ) Z( 1 )) = (^,fci)) USm§l j ' 

Cross product and sum over (^o?^i) give 

[E(Z(°) + Z«) + fc + fci]P(Z j0 = fcb, Zjj = fci | Zj = 1, Zi > 0) 

> minfo, fci}P((Z(°>, Z«) = (fc , fc)). 

This gives the result since P(Zjo = Aft, Zji = A;i | Zj = 1, Z\ > 0) > a with 

= min{fco,fc 1 }P((Z(°),Z( 1 )) = (fco,fc 1 )) 
a E{Z(°) + ZW) + k + k 1 ' 

In the case P(Z^ZW = 0) = 1, we choose (k , h) G N* 2 such that P 2 ((Z , 
^i) = (^0: &i)) > [using (2)]. We make then the same proof as above with 
Zj = 2 and 

min{fco, &i}P2((^o, z i) = ( fc o, h)) 
E 2 (Z + Z 1 ) + fc + A; 1 

so that the result follows as previously. □ 

Thus if P(Z(°)ZW =0) = 1, we need to prove that there are many cells 
with more than two parasites in a contaminated cell line. 



Lemma 4.3. If (3 := P(Z< > > 2 or > 2) > then 

2jo > 2or Zj 



>f P(#{j < i : Zj > 2or Zj! > 2} > /3n/2|Z; > 0) 1. 



Proof. For all i £ G n and j < i, let k such that i = jk, then for every 
a > 0, 

P(Z j0 > 2 or Zj! > 2 | Zj = a, Z K > 0) > P(Z > 2 or Z x > 2 | Z k > 0) > fi. 

Then conditionally on Z\ > 0, #{j < i : Zjo > 2 or Zji > 2} > J2k=o flk, where 
(Pk)i<k<n are i-i-d. and distributed as a Bernoulli(/3). Conclude with the law 
of large numbers. □ 

Proof of Theorem 4.1. We consider first the case when rriQ + mi > 1, 
work conditionally on Ext c and choose i S 5T* . 

If P(Z(°)ZW = 0) / 1, Lemma 4.2 (with k = I) entails that a.s. under 
P(-|Zi>0), 

#{j<i:Z jo >0,Z jl >0} = cx>. 
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Using the branching property and the fact that the probability of nonex- 
tinction of parasites is positive ensures that #G* qq a _ S- 

If F(Z^Z^ = 0) = 1 then P(Z(°) > 2 or ZW > 2) > and by Lemma 
4.3, we have a.s. on P(-|Zj > 0), 

#{j<i:Zjo>2arZji>2} = oo. 



Using as above Lemma 4.2 (with k = 2) and the fact that the probability of 
nonextinction of parasites is positive ensures that #G* r L-3? oo a.s. 

We consider now the case when mo + m\ < 1 and work conditionally on 
i = (ax, . . . , a n ) € G* . We denote ij := (a±, . . . , ay-i, 1 — aj) for 1 < j < n. 

If P(Z(°)Z( 1 ) = 0) / 1, Lemma 4.2 entails that 

(19) VI < j < n, k > 1 P(Z ij . > | %_x = k, Z\ > 0) > a. 

Moreover, if Z\. > 0, then the tree of contaminated cells rooted in ij dies 
out and so has at least one leaf. So by the branching property, the number 
of leaves converges in probability to infinity as n tends to infinity. 

If P(Z(°)ZW = 0) = 1, (19) holds for k > 2 and Lemma 4.3 allows to 
conclude similarly in this case. □ 

5. Proportion of contaminated cells with a given number of parasites. 

We determine here the asymptotics of the number of contaminated cells 
and the proportion of cells with k parasites, defined as 

In that view, we introduce the Banach space Z X (N) and the subset of fre- 
quencies S X (N) which we endow with the norm || • ||i defined by: 



oo 

l-E? I ) 

i=0 ) i=0 



^(N) := < (xi)ieN :E \ Xi \ < 00 f ' IK»i)ieN||i = E 



S X (N) := | (/Oiew : Vi G N, /i € R + , 53 /* = 1 j. 

We shall work conditionally on Ext c or ^ n > and introduce 
(20) P* :=P(- | Ext c ), P n :=P(- | Z n > 0). 

The asymptotics of the proportions depend naturally on the distribution 
of (Z^°\Z^) and we determine five different behaviors according to the 
bivariate value of (mo, mi). 

The proofs of the convergences use the asymptotic distribution of the 
number of parasites of a typical contaminated cell at generation n, which 
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is equal to F n (Zjj n € •), where U n is a uniform random variable in G* in- 
dependent of (^i)igT*. This distribution is different from the distribution 
of Z*, that is the number of parasites of a random cell line conditioned to 
be contaminated at generation n. The following example even proves that 
P n (Z\j n £ •) and F(Z* £ •) could be a priori very different. 

Example. Suppose that generation n (fixed) contains 100 cells with 
1 parasite (and no other contaminated cells) with probability 1/2 and it 
contains 1 cell with 100 parasites with probability 1/2 (and no other con- 
taminated cells). Compare then 



Actually the convergence of (Z*) ne ^ leads to the result obtained by Kim- 
mel [15] in the continuous analogue of this model. That is, 



tends to P(T = k) whereas we are here interested in the expectation of Fk(n). 

A sufficient condition to get the equality of the two distributions is that 
#G* is deterministic, which does not hold here. But in the case when 
(mo, mi) € -D3, we shall prove that #G* is asymptotically proportional to 
(mo + mi) n as n — > 00 (forthcoming Proposition 6.3). This enables us to 
control P n (Z\j n G •) by the distribution of F(Z* £ •). More precisely, it is 
sufficient to prove the separation of descendances of parasites (Proposition 
6.4) and the control of filled-in cells (Lemma 6.5) using the results about the 
BPRE Z*. These two results are the keys for Theorems 5.2, 5.5 and 5.7. Sim- 
ilarly, when (mo, mi) € -D5, we already know that #G* is approximatively 
equivalent to 2™. Then the fact that Z* explodes as n — > 00 (by Proposition 
2.4) will ensure that the proportion of filled-in cells among contaminated 
cells tends to one (Theorem 5.1 below). 

5.1. Case (mo, mi) 6 -D5 (m > 1). In that case, recall that condition- 
ally on Ext c , #G* is asymptotically proportional to 2 n (by Theorem 3.1). 
Moreover the contaminated cells become largely infected, as stated below. 

Theorem 5.1. Conditionally on Ext c , for every k € N, Fk(n) converges 
in probability to as n — > 00, that is, 



P n (Z Un = 1) = 1/2, 
P(Z* = 1) = 100/101 



F n {Z Un = 100) = 1/2; 
P(Z* = 100) = 1/101. 



(21) 



¥(Z n = k) _ E ig G n mz i= k) _ E(#{i £ G n : Zj = k}) 
P(Z n > 0) E ieG „ E(l Zl>0 ) " E(#G* ) 



VK,e>0 



F 



( 



#{i€G w :Zi>if} 



>l-£ 



) 
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If rriQ = mi, the number of parasites in a contaminated cell is of the same 
order as uiq . More precisely, for every e > 0, 

g{p . r ^ r ^ e)Ho . 

Proof. In that case, use Theorem 3.1 and (20) to get that there exists 
a nonnegative random variable L such that 

(22) #G* n > 2 n L, P* (L = 0) = 0. 

Let K, 7} and e > and put B n (K, n) := { #{ie ^ G f i - ir} > r?} n Ext c , then 

E 1 {^i<X}> ? ? 2 "^B n (^r ? ) 
i€G* 

which gives, taking expectations, 

^ nY J ^J- n Hz i <K } ) F(0 <Z n <K) 

K " 1] 7] 

Use then Proposition 2.4 and (22) to choose n large enough so that 

v(K, V ))<e, 



which completes the proof of the theorem. In the case m,Q = mi = m, follow 
the proof above and use that Z n /m n converges to a positive limit on Ext c 
(see [4]) to get the finer result given after the theorem. □ 

5.2. Case (mo, mi) € D3 (m < 1). We assume here E(Z'^ 2 ) < 00 and 
prove that (Fk(n))k£N converges to a deterministic limit. We prove the con- 
vergence thanks to the Cauchy criterion [using completeness of Z 1 (N)]. The 
fact that the limit is deterministic is a consequence of the separation of the 
descendances of parasites and the law of large numbers. Once we know this 
limit is deterministic, we identify it with the Yaglom limit T (see Section 6.1 
for proofs). 

Theorem 5.2. Conditionally on Ext c , asn^-oo, (i 7 fc(^))fceN converges 
in probability in § X (N) to (P(T = k))k^. 

Remark. We get here a realization of the Yaglom distribution T. 

The limit just depends on the one-dimensional distributions of (Z^ , ) . 
More precisely, recall that the probability generating function G of T is char- 
acterized by (15). 

This theorem still holds starting from k parasites. We also easily get a 
similar result in the case when a cell gives birth to N cells [N £ N). 
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As an application, we can obtain numerically the Yaglom quasistationary 
distribution of any BGW. Let Z be the reproduction law of a BGW with 
mean m < 1 and choose N such that Nm > 1. Consider Kimmel's model 
where each cell divides into N daughter cells and = Z^ = ■ ■ ■ = Z^ = 
Z. Computing then the asymptotic of the proportions of contaminated cells 
with k parasites gives the Yaglom quasistationary distribution associated to 
Z. If P(Ext) 7^ 0, one can start from many parasites "to avoid" extinction. 

More generally, we can obtain similarly the Yaglom quasistationary dis- 
tribution of any BPRE with finite number k of environments such that 

This theorem is in the same vein as Theorem 11 in [14]. But we can not 
follow the same approach as Guyon for the proof. Indeed we have to consider 
here the proportions among the contaminated cells in generation n whereas 
Guyon considers proportions among all cells in generation n. Unfortunately, 
the subtree of contaminated cells is itself random and induces long-range 
dependences between cells lines, so that Guyon's arguments do not hold 
here. Moreover, Theorem 11 in [14] relies on an ergodicity hypothesis which 
cannot be circumvented. 

Example. We give two examples when the limit can be calculated. 
Trivial case: P(Z<°) £ {0, 1}, Z« G {0, 1}) = 1 leads to P(T = 1) = 1. 
Symmetric linear fractional case: p £ ]0,1[, 6€]0,(l-p) 2 [ and 

F(Z^ = k)=F{Z (1) = k) = bp k ~ 1 ifjfe>l 

and P(Z(°) = 0) = P(Z^ = 0) = (1 - b - p)/(l -p). Then m = mi = 6/(1 - 
p) 2 < 1 and letting so be the root of fo(s) = s larger than 1, 

W;>1 P(T = k) = (so - 1)/sq- 

As asymptotically we know the number of parasites and the proportion of 
cells with k parasites, we get the number of contaminated cells [recall that 
W is given by (6)]. 

Corollary 5.3. Conditionally on Ext c , the following convergences hold 
in probability 

~Z~ n ' E(T)' (mo + mi)» ~^ E(T) ' 

We can also consider the ancestors at generation n of the cells of dT* , 
which amounts to considering 

Fk{n,p) = 
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and let p — > oo. Letting then n — > oo yields the biased Yaglom quasistationary 
distribution, thanks to the separation of descendances of parasites. 

Corollary 5.4. Conditionally on Ext c , for every k € N, Fk(n,p) con- 
verges in probability in § 1 (N) as p iends to infinity. This limit converges in 
probability in § 1 (N) as n — > oo: 

VfcGN lim lim F k (n,p) = - i~ 

n^oop- >oo J|h,( X J 

We get here an interpretation of the fact that the stationary distribution 
of the Q-process associated to the BPRE (Z n ) ng N is the size-biased Yaglom 
limit (see [1]). 

5.3. Case (mo, mi) € Z?2- In that case, the parasites die out. So we con- 
dition by Z n > 0, we still assume E(Z^ 2 ) < oo and we get a similar result. 



Theorem 5.5. As oo, (Fk(n))k<zfq conditioned by Z n > converges 
in distribution on S X (N) to (P(T = k))^^. 

The proof follows that of the previous theorem. Indeed (13) is still satisfied 
and we can use the same results on the BPRE (Z n ) n ^. There are only two 
differences. First, we work under P n instead of P*. Moreover Z n satisfies now 
¥(Z n > 0) n ~°° 2/(Var(Z(°) + Z^)n) and Z n /n conditioned to be nonzero 
converges in distribution as n — ► oo to an exponential variable £ of parameter 
2/(771 + 1) (see Section 2.1). As above, we can derive the following result. 

Corollary 5.6. As n — > 00, #G* /n conditioned by #G* > converges 
in distribution to £/E(T). 

5.4. Case (mo, mi) £ D\. In this case, the number of contaminated cells 
does not explode and the number of cells of type k at generation n condi- 
tioned by the survival of parasites in this generation converges weakly to a 
nondeterministic limit (see Section 7 for proofs). 

Theorem 5.7. As n — > 00, (#{i G G* : Z\ = fc})fc e ^ conditioned on Z n > 
converges in distribution on / X (N) to a random sequence (N^keN which 
satisfies E(X)feeN kNk) < 00. 

As above, we get: 

Corollary 5.8. #G* conditioned by #G* > converges in distribu- 
tion to a positive finite random variable. 
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Picking a cell uniformly on dT* leads again to the size-biased distribution. 



Corollary 5.9. For every n e N, (#{i G G* n+p :Z-^ n = k})ken condi- 
tioned on Z n+p > converges weakly in Z 1 (N) to a random sequence as p 
tends to infinity. This limit converges weakly as n — > oo . 

kN k 



VA;GN 



lim lim ^{i G 

n-^oo p— >oo 



"n+p ■ ^i\n 



k}\Z n+p >0- 



5.5. Remaining domain: (mo, mi) £ D 4. In this domain, the asymptotic 
of the mean of the number of contaminated cells, that is E(#G*) = 2 n ¥(Z n > 
0), is different from the previous ones. 

Recalling Section 2.2, this asymptotic depends on three subdomains, the 
interior of D4 and the two connex components of its boundary. More pre- 
cisely, it depends on m^mi = 1 or m^mi < 1 and molog(mo) + milog(mi) 
is positive or zero. 

If (mo, mi) G -D4 and mo < 1 < m±, using (17) and a coupling argument 
with Corollary 5.3, one can prove that 



sup<^ P 



net 



2 n F{Z n > 0) 



> A 



(m + m ) 



0. 



where mo = (1 + J 1 + 4(mo — mg))/2 > 1. Thus #G* grows geometrically 
and one can naturally conjecture that #G^ is asymptotically proportional 
to E(#G*) = 2"P(Z„>0). 

Moreover separation of descendances of parasites, control of filled-in cells 
and Corollary 5.4 do not hold in this case. Thus determining the limit be- 
haviors here requires a different approach. 

Finally, note that in the subdomain m^mi = 1 (boundary of D5), (Z*) n ^ 
explodes (see [2]) so the asymptotic proportion of contaminated cells which 
are arbitrarily largely contaminated should be equal to 1 as in Theorem 5.1. 

6. Proofs in the case (mo, mi) £ D3. We assume in this section that 
< 00 (i.e., m < 00) and we start with giving some technical results. 

6.1. Preliminaries. First, note that for all u, v G / 1 (N*), we have 



(23) 



Fill 



Mi 



u — v 



Fill 



+ 



Mli 



Mli 



Mli 



Fill 



<2'- 



v \\i 



Fill 



Moreover by (6), there exist two random variables C and D a.s finite such 
that 

(24) VneN C < - Z \ < D a.s., F*(C = 0) = F*(D = 0) = 

(2m) n 
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and as Clneni^n > 0} = {Vn G N : Z n > 0}, we have 

(25) sup{|P n (A) -P*(A)|}™0. 

A 

We focus now on the BPRE (Z n ) n( zfq. First, by induction and convexity 
of f a , we have for every i G G n (see Section 2.2 for the notation) 

(26) P(Z i >0) = l-/ i (0)<m i . 

Then identities (26) and (14) entail that there exists M > such that 

(27) M< f < Z " >0 »<l. 

Moreover, by Corolary 2.3 in [1], we have 

(28) lim sup{E(Z n t Zn > K | Z n > 0)} = 0. 

Finally, following the proof of Theorem 1.2 in [13] (see [7], Section 2.1 for 
details) ensures that, if (Zn^)nen and (Zn^) n& ^ are two independent BPRE 
distributed as (Z n ) ne pj, we have 

P(ZW > 0, Z® > o) = o(P(Z n > 0)) = o(m n ) (n oo). 

Then, we have 

(29) 2~ n > 0) 2 = o{m n ) (n -» oo). 

iGG„ 

6.2. Estimation of #G* . We prove here that the number of parasites 
which belong to filled-in cells is negligible compared to the total number of 
parasites (see also Lemma 6.5 for a result of the same kind). To prove this 
result, we use its counterpart for BPRE (Z n ) ng N conditioned to be nonzero. 



Lemma 6.1. For every t] > 0, 

su P (p* ( gigg^i^l > v ) \K=a? o. 

neN I V Z n 

PROOF. Let 7? > and write 



Then 



«,):={^|^>.}nExf. 



iGG* 
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Using (24), we have 

ieG* 

so that taking expectations, 



rn 



-E(2-« £ Z { t {Zi>K \ >E(l An{KtV) C) V 



m n E(Z n t {Zn>K} )/r] > E(t An ( K>r]) C). 
Then, by (28), we have 

lim sup{E(t An(Ktri) C)} = 0. 

Then observe that Va > 0, inf P »( j4 )> Q {E(Cl yl )} > 0. So 3K > such that 
VK>K , VneN, 

F*(A n (K,7 ] ))<a, 
which completes the proof. □ 

First, for any e > 0, choose K using the previous lemma such that 
'EieG^z,^} > 1/2 \ =1 _ ¥ *f giggS gi^>^} < 1/2 \ > x _ £/2> 

Adding that conditionally on Ext c , Z n qq a . S) gives the following result. 



Proposition 6.2. Let e > 0, i/tere existe if 6 N suc/i that VW 6 N, 
3no € N smc/i i/iai Vn > hq , 



Zi1 {Zi < k} >n\>1-6. 
VieG* / 



Second, we derive an estimation of By Lemma 6.1, the cells are 

not very contaminated so the number of contaminated cells is asymptoti- 
cally proportional to the number of parasites, which is a Bienayme Galton- 
Watson process. 

PROPOSITION 6.3. For every e > 0, there exist A,B>0 such that for 
every n€N, 

G \e[A,B})>l-e. 



(2m) 
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Proof. First use (24) to get 



#g; < < ^ 



(2m) n ~ (2m) n 
Moreover using again (24), we have 

#g; > Ei £ G* z \\z,<k} > c Ei G G* Zii {z .< K} 



(2m) n ~ K(2m) n ~ K Z n 
and Lemma 6.1 gives the result. □ 

6.3. Separation of the descendances of parasites. Start with two parasites 
and consider the BPRE (Z n ) nG ^. Even when conditioning on the survival of 
their descendance, the descendance of one of them dies out. This ensures that 
two distinct parasites in generation n do not have descendants which belong 
to the same cell in generation n + q if q is large enough. More precisely, we 
define N n (i) as the number of parasites of cell i\n whose descendance is still 
alive in cell i and we prove the following result. 

Proposition 6.4. \/K € N, Ve, r\ > 0, 3q e N such that Vn € N, we have 
mSf f#{^G* n+q :Z i{n <K,N n (i)>2} 



PROOF. Let feN,f|>0 and consider for A > 0, 



>r))<e. 



Then 



+<■/ 



l BSW #{i G G; +g : Z i]n < K, N n (i) > 2} > t EUv)V A(2m) n ^ 
so that taking expectations, 

2_(n+9)E (EieG n+9 1 {2 i | n <^,JV„(i)>2}) 



P(^(r?)) < 



< 



< 



i]Am n+q 

2"" E jg G„ < gj < E ig G g MjVo(i) > 2) 

P(Z n >0)2-«Ei e6 Pjr(M)(i)>2) 



r]Am n+ i 

As we have ( 2 ) ways to choose two parasites among K and they both survive 
along i with probability P(^i > 0) 2 , we have 

IPV(iVo(i) > 2) < (^V(^i>0) 2 . 
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Then 



r qAm q 

Conclude choosing A in agreement with Proposition 6.3 and q with (29). 

□ 

6.4. Control of filled-in cells. Here we prove that filled-in cells have 
asymptotically no impact on the proportions of cells with a given number 
of parasites. 

Lemma 6.5. Ve, r] > 0, 3K e N such that Vn, q G N, we have 
m f#{ieG* n+q :Z iln >K} 



Proof. Let rj > 0, A > and consider 



>f? <£. 



then 



- J l(2m)"+9 



l F « ( „)#{i € G* : Z, n >K}> l^ M T,A(2m) n+ ". 



I 1 fc ^n+g • ^i|n -"■ J — A F„(r]) 

Taking expectations leads to 

2-(^)E(E i6Gn+g l { ^ 1|n >ic^ i >0}) 



P(i™fa)) < 



< 



2- (n+9) Ei eGn+ ,P(^i|n> ^,^i>0) 



< 



E fc >/r 2~ n E je G„ P(3 = fe)2-« EieG, Pfc(^ > 0) 



r]Am n +i 

Moreover, ¥ k (Zi > 0) = 1 - (1 - P(Z| > 0)) fc < &P(Zi > 0) and we have 



r)Am n +i 



E(Z n I {M} ) . 

< — *- using (27). 

r\Amr 

By (29), we get 

lim sup{P(F 1 ?( ?? ))} = 0. 
Complete the proof choosing A in agreement with Proposition 6.3. □ 
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6.5. Proof of Theorem 5.2. Consider the contaminated cells in gener- 
ation n + q. Their ancestors in generation n are cells which are not very 
contaminated (by Lemma 6.5). Then taking q large, the parasites of a con- 
taminated cell in generation n + q come from a same parasite in generation 
n (separation of the descendances of parasites, Proposition 6.4). Thus at 
generation n + q, everything occurs as if all parasites from generation n be- 
longed to different cells. As the number of parasites at generation n tends to 
infinity (n — > oo, mo + nii > 1), we have a law of large numbers phenomenon 
and get a deterministic limit. 

Step 1. We prove that for all e,r] > 0, there exist no € N and / 6 S 1 (N) 
such that for every n > no, 

r , {\\{F k {n)) k&i -f\\x>Tj)<e. 

For every k G N* and every parasite p in generation n, we denote by Y k (p) 
the number of cells in generation n + q which contain at least k parasites, 
exactly k of which have p as an ancestor. By convention, Y$(p) = 0. That 
is, writing for p parasite, p <^-> i when p belongs to the cell i and p|n its 
ancestor (parasite) in generation n, 

Y k(P)= E 1 #{r:r^,r|n=p}=fe ) & € N*. 

By the branching property, (Y k (p))keN [p £ ^( n )] are i-i-d. and we denote 
by (5j?)fceN a random variable with this common distribution. Denoting by 
Vk{ji) the set of parasites in generation n which belong to a cell containing 
at most K parasites, we have 



#{iGG; +g :Zi = fc}- Y, n 9 (p: 

peV K (n) 



E 

(30) 

< (K + l)#{i G G* n+q : Z i]n < K, N n (i) > 2} + #{i G G; +g : Z i{n > K}. 
Indeed, the left-hand side of (30) is less than 

+ #{ieG* n+q :Z iln >K}. 



E 

fceN* 



#{\eG* n+q :Z i = k,Z. An <K}- J2 Y kiP) 

pePx(n) 

And recalling that N n (i) is the number of parasites of cell i|n whose descen 
dance is still alive in cell i, we get the following equalities: 

E Y k(P)= E E 1 #{r:r->i,r|n=p}=fe 

and 

1 Z i =k,Z iln <K,N n (i)=l = 1 JV„(i)=l E 1 #{r:r^i,r|n.=p}=fc 

pev K (n) 
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which ensure 



E 

fceN* 



#{ieG; +g :Z i = A;,Z i | ri <K}- Y k& 

P ev K (n) 



<E E 

fcGN* ieG n+9 ,Af„(i)>2 



^-Zi=k,Z i]n <K ~ E 1 #{r:r<-i,r|n=p}=fc 
peP/c(n) 



< #{i € G* +9 :Zi| n < K, N n (i) > 2} + l#{r : r^i,r|n=p}>0 

ieG„ +9 ,7V„(i)>2 

<#{i6G; +? :Zi| n <ir, jV n (i)>2}+ Kl z i]n <K 

ieG n+9 ,7V„(i)>2 
= (K + l)#{i G G; +g : Z- An < K, N n (i) > 2}. 

We shall now prove that the quantities on the right-hand side of (30) are 
small when n and q are large enough and that J2peT K M Y k (p) f° now a 1& W 
of large number. To that purpose, let s, r] > and for all K,k,n,q>0 define 



Spg-p K (n)^fc (p) 

SfceN J2 p eV K (n) Yk (p) 



First, by Proposition 6.2 and (25), G N such that ViV G N, 3n x G N 
such that MK > K±, Vn > rii, 

(31) V n {\V K {n)\>N)>l-e. 
Moreover by Lemma 6.5, 3K 2 > -fTi such that \/n,q G N, 

(32) r ^ ^_^n 2/ 

And by Proposition 6.4, 3qo G N such that Vn G N 



>?? <£• 



/#{iGG; +w :Z i |„<^ 2 ,iV„(i)>2} 



>r]/(K 2 + l)) <e. 



(33) 

\ ft^n+qo 

Use then (30), (32) and (33) to get 

/EfceN* € G n+(?0 : Z ; = k} - T, P er K2 (n) Y k°(P 

Then by (23), for every n G N, we have 

(34) F*(\\(F k (n + g )) fceN - (Cf (n, gb)) fceN |li > 4iy) < 2e. 



> 2rj < 2e. 
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Second, conditionally on Z n > 0, Y k qo (p) [p G 7V 2 ( n )] are i.i.d. Then the 
law of large numbers (LLN) ensures that Vk £ N, as n and so Vki ( n ) becomes 
large: 

G% 2 (n, g ) — ► fk(qo) where f k (q ) := — E ^ ^ 

To see that, divide the numerator and denominator of G^ 2 (n, qo) by #Vk 2 ( n )- 
More precisely, by the LLN, there exists N > such that for all n € N, 

P"(||(Gf 2 (n,g )) fce N* - /Mill >V,V K2 (n) >N) <e. 
So using (31), there exists n\ € N such that for every Vn > m, 

P"(||(Gf 2 (n,g W -/(?o)||i>»7)<2e. 
Finally by (25), there exists 712 > ni such that for every n > rt2, 

(35) P*(||(Gf 2 (n,g )W -/M||i >t?)<3 £ . 

As a conclusion, using (34) and (35), we have proved that for all £,77 > 0, 
and for every n>ri2 + qo, 

r(||(F fc (fi)) fceN .-/( go )||i> 577)<3e. 

5*tep 2. Existence of the limit. 

For every I € N, there exist no(l) E N and f(l) £ S 1 (N) such that for every 
n > no(0 

P(||F(n)-/(/)|| 1 >l/2 m )<l/2'. 

Then for all I, V such that 2<1<V: \\f(V) - f{l)\\ x < 1/2' and completeness 
of Z X (N) ensures that (/(O)zeN converges in S 1 (N) to a limit /. Moreover, 
11/(0 — /111 — 1/2' so for every n > no(Z), 

P(||F(n) - f\\ x > 1/2') < 1/2' 

which ensures the convergence in probability of (Ffc(n)) ng N to / as n — > 00. 
.Step 3. Characterization of the limit as /& = P(T = /c). 
By Proposition 2.2, we have 

(36) VfceN P(Z n = jfe I Z n ^ 0) n -=^ P(T = jfe). 
Moreover, for every fcsN*, using (21), 

E(#{i € G n : Zi = k}) E(F k (n)#G* n ) 



¥{Z n = k\Z n ^0) 



E(#G*) E(#G*) 
As Fk(n) converges in probability to a deterministic limit fy, we get 
(37) V7c G N P(Z n = fc I Z n ? 0) ™ / fc . 
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Indeed, by Proposition 6.3, there exists A > such that 

E(#G*) 



(2m) r 



> A 



Then for every rj > 0, using \Ff.(n) — f k \ < 1, we have 



E(F k (n)#G n ) 



E(#G*) 



/a 



< 



E(#G^ I F fc (n) - / fc | l{|F fc (n)-/ fc |Q7}) 

E(#G*) 

E (# G n 1 {|F fc (n)-/ fc |>rj}) 



+ 



<7] + 



E(#G*) 

E (- 2: nl{|F fc (n)-/ fc |>» ? }) 

^(2m) n 



By (11), Z n /(2m) n is bounded in L 2 and it is uniformly integrable. Then, 
thanks to the previous steps, the second term in the last displayed equation 
vanishes as n grows and we get (37). Putting (36) and (37) together proves 
that / fc =P(T = fc). 

6.6. Proof of corollaries. 

Proof of Corollary 5.3. Recall that E(T) < oo (Proposition 2.2) 
and note also that for every K £ N* , 



Then using EieG* Z \^{z-,<K} < z n gives 



E(T) 



1 



< 



1 



E(T) 



1 



EtikF k (n) E(T) 



+ 



1 


EigG* Z i^{Zi<K} 


E(T) 


Zn 



Let r/,e > 0. We use Lemma 6.1 to choose if GN* such that 

'EiGG* Zil{ Zi <K} 



VneN 



1 



>l-7?)>l-£; 
1 



< rj. 



E(T1 T < X ) E(T) 
Choose no G N using Theorem 5.2 so that for every n > hq, 

1 1 



Ek=ikF k (n) E(T1 T < A ') 



<V > 1 
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Then for every n > hq , 



E(T) 



>2 v+¥ L_ v] <2e, 



which proves the convergence in probability of ^G*/i? n to 1/E(T). The 
second convergence follows from (6). □ 

Proof of Corollary 5.4. We write for n,p, k G N, 

#{ieG* n+p :Z iln = k} = (2m)P ^ #{i€G* +p :i|n=j} 

Conditionally on Zj = k, by Corollary 5.3 and separation of descendances of 
parasites, we have the following convergence in probability 

#{i€G; +p :i|n=j} p-«o 
(2^ > 

where Wfc(j) is the sum of k i.i.d. variables distributed as W/E(T). Then, 
using also (6), 

(38, nwm- mw) - k 



E(T) E(T) 

Using again Corollary 5.3, we get the first limit of the corollary 

#{i£G; +r % = fc} p E(T)Ej e G-zj=fcWfcCi) 

A™ #G^ +p W (2m)- 

Moreover, Theorem 5.2 ensures that 

#{j€<G£:Zj = fc} Z ra ^ 

(2m)« fcW (2m)« ~^ E(T) Jfe ' 

And conditionally on #G* > 0, (j G G*) is i.i.d. by the branching 

property and #G* tends to infinity. So the law of large numbers and (38) 
ensure that 

E(T)EjgG*:Z J =fc^fc(j) 

n^o W (2m) n 

E(T) #{j G G* :Z- } = k} E jeG * : z- 3 =k W k (j) p. fc/ fc 



lim 



n^oo (2m)" #{j G G* :Zj = k} E(T) : 

which ends the proof. □ 
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7. Proofs in the case (mo, mi) G D\. We still assume < oo, 

the proof is in the same vein as the proof in the previous section and use the 
separation of the descendances of the parasites. The main difference is that 
Z n does not explode so the limit is not deterministic and the convergence 
holds in distribution. 



Lemma 7.1. For every K > 0, there exists qo G N such that for all q > qo 
and n G N, 

P" +9 ({i G Gl +q : N n (i) >2}^0,Z n < K) < e. 

Proof. Denoting by the event 

{{i G G* n+q : N n (i) >2}^0,Z n < K}, 

we have 

1 E* ^ Yl 1 {N n (i)>2,Z n <K}- 

Thus we can follow the proof of Lemma 6.4. 



DTL+q 



™< V F(N n (i)>2,Z n <K) 



ieG„ 



+1 



¥{Z n+q > 0) 



E ieGn+g P(iV n (i)>2,Z i|n <if) . 

* using (7) 

< P(0 < Z w < K)2-« Ei gGq MjVi) > 2) 
Um n +i 

< 7J^ (27) • 

Conclude with (29). □ 

Proof of Theorem 5.7. 

Step 1. We recall that V n is the set of parasites in generation n, follow 
Step 1 in the proof of Theorem 5.2 and use its notation. Thus, we begin with 
proving that for every e > 0, there exists no G N such that for every n > no, 



w+qi 



|(#{i G G* n+q : = k}) km - (N k (n, q))^ + 0) < e, 



where for all n,q,k>0, N k (n,q) :=J2 P eV(n) Y k(p)- 

First, by (10), there exist K, go G N such that for every q > qo, 

(39) lim F n+q (Z n > K) < e. 
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By Lemma 7.1, there exists q\ > qo such that for every n G N, we have 

(40) P n+<?1 ({i G G ? ; +gi : iV n (i) > 2} / 0, Z n < K) < e. 
And by (39), there exists no > such that for every n > no, 

P n +^(Z n >K) <2e. 

Then 

P" +yi (#{i G G; +91 : iV„(i) > 2} + 0) < 3e. 

Moreover, 

#{iGG; +?1 :jV re (i)>2} = 

=> (#{i G G* n+qi : Zi = k}) km = (N k (n, qi)) ke n- 

Then for every n > no , 

lT +9l (ll(#{i G G* n+qi : Z; = k}) km - (JVfc(n,gi)) fc6 N||i ^ 0) < 3e. 

S'tep 2. As / 1 (N) is separable, we can consider the distance d associated 
with the weak convergence of probabilities on / 1 (N). It is defined for any Pi 
and P2 probabilities by (see Theorem 6.2, Chapter II in [18]) 

d(P 1 ,P 2 ) = sup| J fiw^idw)- J /(«7)F 2 (d«;)| : H/lloo < 1, II/IIlipb < 1 
where 

ll/l| LiP s = sup( ll/( | f- / | ^ l|l :x,,G^(N ) ,^,|. 

We prove now that for every I > 1, there exist no(Z) G N and a measure //(/) 
on N* such that for every n > no(l), 

(41) d(P n ((#{i G G* : Z K = k}) keN G < ■ 

For that purpose, let I G N. By Step 1, choose q, uq G N such that 
Vn > n d(P n+9 ((#{i G G* n+q : Z K = k}) km G •), 

(42) 

F n+ i((N k (n,q)) keN e-))<l/2 l+1 . 

Recall that (Y fc 9 (p)) fcGN (p G V(n)) is an i.i.d. sequence distributed as (Y^ )fceN 
and #V(n) = Z n . Thus, under p n +9 5 N k (n,q) is the sum of Z n variables 
which are i.i.d., distributed as Y k q and independent of Z n , conditionally on 

EfcGNEpg^(n) Y k (p) > 0- 

Moreover ¥ n+q (Z n G •) converges weakly as n — > 00 to a probability v 
[see (9)] and we denote by J\f a random variable with distribution v and 
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by (Y^(i))k£fq(i € N) an i.i.d. sequence independent of TV and distributed as 
(Y^fceN- Then we have for n large enough, 

(43) d(F n +i((N k (n,q)) ke N 6 -),fi(l)) < 1/2', 

where is the distribution of (Za<i<A/'^fe 9 (0)fceN conditionally on 

£fc£NEi<i<A^« > 0- Combining (42) and (43) gives (41). 

Conclusion. As I (N) is complete, the space of probabilities on Z X (N) en- 
dowed with ti is complete (see Theorem 6.5, Chapter II in [18]), (//(Z))z e pj 
converges and we get the convergence of Theorem 5.7. 

We now prove that E(^ fcgN * kN k ) < oo. For all n,K > 0, we have 



mi) 



¥(Z n > 0)K 



E( 2 k#{ieG* n :Z i =k}\Z n >o\ < E(iJ n l { ^>X}!^n > 0) < 

which converges uniformly to as K — > oo using (11). Moreover, Theorem 
5.7 and £;#{i € G* : = fe} < Z n ensure that 

n limE( J2 k#{ieG* n :Zi = k}\Z n >o\ =E[ £ fcJV fc ) . 
rwoo y 1<fc<x y \i<fc<_ft: / 



(^fc#{iGG;:Z i = fc}|^ ri >0j ^E( ^ 

VfceN / VfceN* / 



Thus we get the expected limit 
E 

and recalling Section 2.1, we have also 

E I J2 G G n : Z i = k } \Zn>0)= E(Z n \Z n >0) ™ B'{1) < oo. 

VfceN* / 

This completes the proof. □ 

The proofs of the corollaries follow those of the previous section. 
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